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Abstract. I discuss in parallel two universal phenomena: the independence of statistics of the 
heat capacity and entropy of ideal gases of the same, constant, density of states, on one hand, 
and the independence of statistics of the heat and entropy transport through one-dimensional 
channels, on the other hand. I show that there is a close similarity between the microscopic 
explanations of each of these phenomena. 



1. Introduction 

The thermodynamic equivalence is a name coined by M. H. Lee [1] for the property of ideal Bose 
and Fermi gases confined in two-dimensional (2D) boxes to have the same heat capacity, Cy, 
and entropy, S. This implies that from the thermodynamic perspective, 2D Bose and Fermi 
gases under canonical conditions (i.e. has a constant volume and can exchange only heat with 
the environment) are identical. 

This property was first proved by Auluck and Kothari in Ref. [2] and then rediscovered by 
May [3] and Viefers, Ravndal and Haugset [3]. Nevertheless, it did not receive much attention 
until Lee introduced a unitary description of ideal Bose and Fermi gases [H El [6l [7| in terms of 
polylogarithmic functions [H [9] . 

In Ref. [To], Haldane introduced an extension of the Pauli exclusion principle, which was 
latter called the fractional exclusion statistics (FES). In FES, the system is divided into species 
of particles, that we shall identify by an index, e.g. i = 0,1, . . .. In each species there are 
Gi available single-particle states and Ni particles. The FES is defined by the fact that if we 
change for example Ni by 6Ni, then the dimension of species j, for any j = 0,1, . . ., changes by 
5Gj = —oijidNi, where aij are called the FES parameters. The diagonal parameters [an) are 
called direct parameters, whereas the rest of them [uij , i ^ j) are called the mutual parameters. 
For example for a continuous (i.e. macroscopic) FES system the species of particles are defined 
by coarse-graining the phase space, each grain representing a species of particles [HI [121 fT3] . 

The general properties of the FES parameters were deduced only relatively recently 
[m [151 [E] • A quite general ansatz for the FES parameters that obey the properties introduced 
in [H] is [m [HI 

(^ij = a\f +a\'^6ij. (1) 
The parameters a^j are "extensive", i.e. they are proportional to Gi - the dimension of the 



aij' = aijd, (2) 



species on which they act, 

r. 

(s) 

whereas the parameters a] are direct parameters and are not extensive. 

Typically, in the literature we find exclusion statistics parameters of the (s) type (see e.g. 
Ha da Eni EH Ea Ea El ESD, so in general Oij = for any i ^ j. If, moreover, we impose 
an = a for any i, the thermodynamic calculations simplify considerably. If a = we obtain the 
Bose statistics, whereas if a = 1 we obtain the Fermi statistics. 

If aij = a5ij, one obtains a very surprising result: in 2D systems Cy and S are independent 
of a and equal to the heat capacity and entropy of ideal Bose and Fermi gases. Therefore 2D 
systems of any statistics, i.e. Bose, Fermi or FES, are thermodynamically equivalent |24j . 

The property which leads to the thermodynamic equivalence of 2D ideal gases is the constant 
density of single-particle states (DOS). A system of ideal particles in a d-dimensional (dD) box 
has a DOS of the form 

a,ie) = V ( SV^' - CeW'^-\ (3) 



(27r)'^/2r + 1) Uv 



where e is the single-particle energy, V si the dD volume, m is the mass of the particle and F is the 
gamma function. If d = 2, then o"d=2(e) = Vm/ (vr^^) = C is a constant. We shall see in sectionEl 
that the ideal systems of any statistics, of the same, constant DOS are thermodynamically 
equivalent. 

The thermodynamic equivalence has a correspondent in the transport properties of ideal 
gases. At first, Rego and Kirczenow observed that in the low temperature limit, the heat 
conductivity of ideal particles confined in a ID channel is independent of the dispersion relation 
[26] and, moreover, is independent also of their statistics [27]. The statistics independence of the 
heat conductivity implies the independence of statistics of the entropy current in ID channels, 
as it was shown by Blencowe [28] . 

The independence of the low temperature heat conductivity of the dispersion relation have 
been confirmed experimentally in ID phonons [5S] and photons [3^ [51] channels. 

The low temperature results of Rego and Kirczenow have been extended in Ref. [35] . In this 
paper we pointed to another class of universality in systems of ideal particles, namely we showed 
that the low temperature statistics independence of the ID channel heat conductivity can be 
extended to any temperature. In other words, the heat conductivity of ideal particles through 
a ID channel is independent of the statistics of the particles at any temperature. 

In principle all the universal physical properties should have also an universal (and hopefully 
very simple) microscopic physical explanation; there should be a key, basic feature, common to 
all these systems, that is leading always to similar results. For the thermodynamic equivalence 
the key feature is the spectrum of excitations, which is independent of statistics [23], whereas 
for the ID universality the key feature is the current of excitations, which is also independent of 
statistics [32]. The two microscopic explanations are very closely related and, as we shall see in 
the next two sections, the quantities related to the ID transport have mathematical expressions 
very similar to the ones that describe the equilibrium thermodynamics. 

In section [2] we shall discuss the thermodynamics of systems of constant DOS and we shall 
present the main results and formulae. We shall also present the microscopic explanation of the 
thermodynamic equivalence. 

In section[3]we shall present the calculations of the heat, particle and entropy transport in ID 
channels and, by employing the formulae of section [2l we shall show that they are independent 
of statistics at any temperature. The microscopic interpretation of this phenomenon is given in 
the end of section [3] and is very similar to the microscopic interpretation of the thermodynamic 
equivalence. 



Section [His reserved for conclusions. 



2. Thermodynamic equivalence 

Let us assume that we have a macroscopic system of ideal FES particles, with a DOS of the form 
cr(e) = Ce^. The FES parameters are aij = aSij, for any i and j. We associate the same energy 
to all the particles in a species (e^, i = 0,1, . . .) and the same chemical potential, /i, to all the 
species. If we denote by Gi the number of states in the species i when there are no particles in it, 
then in the presence of Ni particles the number of available states becomes Gi — aNi [T8l [2T] . 
The number of particle arrangements that we can have in the species i is 

Therefore the total number of microscopic configurations in the system, when the population of 
each species is fixed, is given by the product [18] 



The grandcanonical partition function of the system, at temperature T and chemical potential 
/u, is 



(6) 



Z=Y. Z{{G„Ni}) = W{{G„N,})exp ^/37Vi(/i 

{N.} {TV.} L i 

where /3 = l/{kBT) and Z({Gi, Ni}) is the "partition function" corresponding to a particular 
choice of the set {Ni}. 

The standard procedure to find the particle populations of the species, denoted by Ui = 
Ni/Gi, is to maximize Z{{Gi, Ni}) = Z{{Gi,ni}) with respect to the n^'s: dZ{{Gi,ni})/dni = 
0. In this way it was obtained 



n(e) = {w{Ce) + a} \ (7a) 

where w is a function that satisfies the equation 

u;(Ce)"[l + w{Q]^-" = = e'^(^-'^) (7b) 

Now we apply the results above to our system of DOS cr(e) = Ce*. From ^ we write the 
expression for the grandcanonical thermodynamic potential, 

- PV = n = -kBTlog{Z) = knT de C7e' log {(1 - an)/[l + (1 - Q)n]}. (8a) 

Jo 

Integrating by parts and using the equations ([7]) we obtain 

1 r°° JJ 
PV = —- de Ce'+'n{e) ee — - , (8b) 
s + 1 Jo s + 1 



where U is the internal energy of the system. 
Similarly, the total number of particles is 



/ deGe'n{e). (9) 
Jo 




Figure 1. The statistics independent functions, yo (a) and — ep) (b), vs the scaled 
temperature, k-oTa/N . 



All these functions may be calculated for example by expressing e in terms of w, 

/3e = log[u;°(l + u;)^~"] +/3/U, (10) 

but the integrals cannot be performed analytically for general s and T. 

For s = (o" = C), all the thermodynamic quantities can be expressed in terms of elementary 
or polylogarithmic functions [21] . I start with N , which, by simple algebra may be written as 

(11) 



N = kBT a log {I + yo), 
where yo = l/i(^(Ce=o) satisfies the equation 

(i + yo)'^7yo = C=V 



(12) 



We observe from eq. ()lip that yo is a function of the dimensionless variable x = k-eTa/N, 
and does not depend on a. Moreover, as seen in figure [T^, yo(x) is a monotonically decreasing 
function, with yo(x — )■ oo) = and yo(x — )■ 0) — )• oo. 

From the equations pT|) and follows 



exp [(^ - aN/a)/kBT] = 1 - exp [-N/{akBT)], 



(13) 



where we can identify the (generalized) Fermi energy as ep = liuiT^o /x = aN/ a and observe 
that // — ep is also independent of a (see figure [Ija). 
After some algebra we obtain 



n = u 



-{k^Tfa 
{kBTf<y 



1 — a 



log7l + yo) + Li2(-yo) 



a 



log2 (1 + yo) + Li2 



yo 



1 + yo 



(14) 



where Li2{x) = J2T=i^'^/^'^ is Euler's dilogarithm [8]. The two expressions for U in are 
connected by Landen's relation, Li2{x) + Li2[—x/{l — x)] = — (l/2)log^ (1 — x), which is valid 
for any x < 1 [S]. 

Using pip into p4p we obtain 



-(^f^)^.(-»). 



a 



yo 
i + yo 



1 

2" ~ 2 

a 

+ - 



(15a) 
(15b) 
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Figure 2. The specific heat (a), cv/Z^b = Cy /{NkB), and entropy per particle (b), s//cb = 
S/{NkB), for a system of particles of constant DOS, a, and any a. 



In the expression ()15ap . Up = —{kBT)'^aLi2{—yo) is the internal energy of a Fermi gas {U = Up, 
if a = 1), whereas in the expression (|15bp . Ub = {kBT)'^crLi2[yo/{l + yo)], represents the internal 
energy of a Bose gas {U = C/b> if « = 0). 

Since yo does not depend on a, but just on k^Ta/N and T, it is obvious that the differences 
between the thermodynamic potentials of gases with different a's come just from an additive 
constant. All the temperature dependence is the same. 

To calculate Cy, we use the expression 

and obtain 

Cv = -^^^ - IklTaLi^i-yo). (17a) 

To calculate the entropy we use S = \J jT + PV/T — ^N/T, together with equations (fTT]l . 
([12]) and ([E]), to obtain 

S = -klTa[2Li2{-yo)+log{l + yo)logyo]. (17b) 

Equations (jl7ap and (|17b|) are independent of a and they express eventually in the most compact 
way the thermodynamic equivalence. 

In figure [2] (a) and (b) we plotted cy/ZcB = Cy /(A^/cb) and s/k^ = S/{NkB), respectively. 

Another way to calculate the entropy is by using the definition 

TS = kBTlogW, (18) 

where W is given by ([5]). From the expression ([5]) we take only the most probable distribution, 
i.e. the one that satisfies the equations ([7]), and, after some algebra we obtain 

logW = {[1 + (1 - a)ni] log[l + (1 - a)ni] - njogni - [1 - arn] log[l - arii]} 

i 

= a de{[l + (1 - a)n(e)] log[l + (1 - a)n{e)] - n{e) log n(e) 
Jo 

-[1 - an(e)]log[l - an(e)]} (19) 
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Figure 3. Excitations in a Bose and a Fermi gas, each of 10 particles and the same DOS. To 
each configuration of bosons it corresponds a configuration of fermions, with the same excitation 
energy. In (a) the two systems are in the ground state; in (b) six particles in each are excited 
on the first "level"; from the first level, 3 particles are excited two levels up (c), whereas in (d), 
one particle from the uppermost "level" is excited two states up. 



where we have used the Stirling approximation, log(A^!) ~ A^log(A^/e), valid for ^ 1, and 
neglected the terms of the order of . Plugging ([7]) into (fT9|) we obtain 



TS 



[w + 1) \og{w + 1) — It? log If; 



JO 



oo 

yo 



w + a 
log{w + 1) log w 



de 



w 

logjy + 1) 

y 



w + l 
logy 



de 



de 



-(fcBT)V[2Lz2(-2/o) + log (1 + yo) logyo], 



(20) 



which is equation ()17bp . as expected. 

Since according to Eq. (|lip yo ^ oo when T — ?• 0, making use of the asymptotic behaviour 
of the dilogarithm, Li2{—yo) ~ — [vr^/6 + log^(yo)/2] one can show that in the limit of low 
temperatures 

Cy - (7rV3)A;|ra{l - 0[log2(l + yo)/yo]}. (21) 



2.1. Microscopic interpretation 

The microscopic interpretation of the remarkable identities expressed by Eqs. (jl7p was given 
in Ref. [23]. We showed in that paper that if we have two FES systems of the same, constant, 
DOS, but different a's, one can establish a one-to-one correspondence between configurations 
of particles in the two systems that have the same excitation energy. The excitation energy is 
the difference between the energy of the system in the given configuration and the energy of the 
ground state; the average of the excitation energy is Ub, defined above. 

In figure [3] we give an example of correspondence between configurations of fermions and 
configurations of bosons with the same excitation energies. In this figure we can see how the 
correspondence is established. Based on this method, in Ref. [33] we showed how we can 
transform a general system of fermions into a system of bosons. 
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Figure 4. Schematic view of the heat and particle transport experiment. The system, which 
consists of two reservoirs connected by a quasi-lD wire, contains FES particles of parameter 
a. The two reservoirs are at the temperatures and chemical potentials (Ti,/ii) and {T2,ij,2), 
respectively. 



3. Transport equivalence 

In Ref. [27], Rego and Kirczenow proved theoretically that the low temperature limit of the 
ID heat conductance if ideal particles is equal to kq = Tr'^k^^T / {?>h) and is independent of 
the statistics of the particles. We extended this result in Ref. [32] and showed that the heat 
conductivity of ID channels is independent of statistics at any temperature. 

Let us now consider a two-terminal transport experiment in which the reservoirs 1 and 2 are 
connected by a quasi ID wire, like in figure [H The temperatures and chemical potentials in the 
two reservoirs will be denoted by and /x^j respectively (r = 1,2). The particles in the system 
obey FES of parameter a. If there are Mc transport channels through the wire, the Landauer 
formula for the particle and energy fluxes between the two reservoirs are \n\ [28] 



^ = E/ ^Vp{k)[ni{k)-n2{k)]Uk), (22a) 
^ = T.j :^ev,mni{k)-n2{k)]Uk). (22b) 



respectively, where Vp{k) is the group velocity of the particles of momentum k and nr{k) denotes 
the thermal particle population in the reservoir r, i.e. nr{k) = n\er{k) ^ jir ,Tr\ and satisfies Eqs. 
d?]); finally, C,p{k) is the particle transmission coefficient through the wire. 

Following [32] and using the thermodynamic relation for heat exchange, Q = T6S = 5U—fi6N 
(see in section [2] the relation for S), where 6S is the entropy variation in the process, we write 
the heat flux, 

Q = U -III (23) 

For massless particles, like phonons, U = Q, since /x = 0. 

Since the particle group velocity is Vp{k) = h^^ {de{k) / dk) , and assuming further that 
Cp{k) = Cp is independent of k, Eqs. (p2]) get the simple form 

I = Et/ de[n{e,fii,Ti)-n{e,fX2,T2)] (24a) 



^ = Et/ d«N(e,m,ri)-^e,M2,T2)] (24b) 

where ep(0) is the lowest energy level in the channel p. 

We are working in the linear regime, i.e. |Ti — T2\lTr <C 1 and — ^2! ^ ^B^^r (^ = 1) 2), 
so we shall denote the average temperature, (Ti + T2)/2, by T. We introduce the notations 
h{f^,T) = "EpIpAf^iT) and ifi{fi,T) = J2pUp,i{n,T), with 



Ip,iif^,T) = ^ naie,ti,T)d€, (25a) 
Up4fi,T) = ^ / en«(e,/i,r)de, (25b) 



and define / and U as 



/ - ^^^ = ^^T + ^5f^, (26a) 

Rego and Kirczenow calculated the low temperature limit of n^y = dUi/dT - which we shall 

call here the energy conductance - and obtained ~ A/"cKo for all a's [STIES], if ep(0) takes 
the same value for all p. 

If we use the relations ([23|) and (p6]) . together with the condition / = 0, we obtain [32] Q = U 
and 

We can observe now already that in general the expression (j27p for the heat conductivity is 
similar to the expression for the heat capacity, if we replace Ui by U, i.e. the energy current by 
the internal energy. 

Let us now set M = 1 and focus only on one channel conduction. In such a case, 
Ii(/U,r) = Ip=i^i{n,T) and we observe that if we set ep=i(0) = 0, then Ii(/x,T) and Ui{fj,,T) have 
the same expressions as the particle number, N{fj,,T), and internal energy, U{fi,T), respectively, 
of a FES system of parameter a and constant DOS, a = Cp=i/h, in equilibrium at temperature 
T and chemical potential fi. Therefore we can apply here directly the results from the previous 
section. 

Let us drop the subscript 1 from Ii{fi,T) and Ui{fi,T) and rewrite the equilibrium 
thermodynamics equations of section [21 using the quantities corresponding to the stationary 
transport. We define again a quantity i/q, which satisfies 

/(/.,T) = ^log(l+yo), (28a) 
(l + 2/o)'~7yo = e-'5^ (28b) 

exp [(// - ahI/C)/kBT] = 1 - exp [-hi /{(kBT)], (28c) 
in analogy to eqs. (fTT]) - (fT3]) 



Prom the eqs. ()28p and observing that U{fi, T) has an expression similar to that of C/(m, T), 
we obtain 



2 



2 C h 2C h Vyo + 1/ 



where we used again the Landen's relations. 
From Eqs. ()28p and using 



and 



we calculate 



dyo\ _ fj-yojl + yo) , . 

dTj^ kBT^{l + ayo) ^^^^ 



kBT{l + ayo) 



Ofj, h[l + ayo) 



dU 
dT 
dU_ 

djjL 



-(1 - a)^^ - 2klT^Li2{-y^) + 



hi dl 
hi dl 



-(1 - a)^^ - 2klT^Li2{-yo) + 



h 

c 
h' 



^^(1 + yo) 
r(l + ayo)' 

1(1 + yo) 
1 + ayo 



(31a) 
(31b) 



With eqs. ([30]) and ([3T]) and using the identity 







a/ 




^ 91 






7=0 


~ dT 
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J 



(32) 



we can calculate the heat conductivity. Plugging (|32|) into (|27j) we obtain 

-1 



at/ af/ ai 

df~~d^Idf 



a/ 
a^u 



(33) 



which, if we replace I hy N and U by C/, becomes identical to the expression for the heat capacity, 
Cv (|17ap . Therefore we can transcribe Eq. ()17ap into an equation for the heat conductivity as 



I^h 1 + yo 2klTC 
TC yo h 



Li2{-yo), 



(34) 



which is independent of a at any temperature. In the low temperature limit Eq. p4 
becomes the universal asymptotic expression, kq, which is equivalent to ()2ip from equilibrium 
thermodynamics, but with o" = C/^- 

The heat flux is related to the entropy flux by the relation dUij=o = TdSi, from where we 
obtain an entropy flux 



S 







which is consistent with the definition used in equilibrium thermodynamics, 
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Figure 5. The microscopic analysis of the statistics independence of the heat conductivity: there 
is a one-to-one correspondence between micro-configurations of particles of different statistics 
which have the same excitation energy, UB,n,i, and therefore carry the same heat fluxes. 

But since S of (|35bp is consistent with the Boltzmann definition (jlSp and therefore with the 
integral of eq. ()20p . it follows that the definition ()35ap is also consistent with 

S = ^J^^{[l + {l-a)n{e)]log[l + {l-a)n{e)] 

— n(e) log[n(e)] — [1 — an(e)] log[l — an{e)]}de 
k'^TC 

= -^^[2Li2(-?/o) + log(l + yo)log2/o], (36) 

which represents the fiux of the number of configurations of particle populations, {n(e)}, from 
one reservoir to the other. But we notice again that eq. ([36]) . like eq. ([201) . is independent of a 
at any temperature. Therefore the entropy fiux is independent of a at any T. 



3.1. Microscopic interpretation 

From the r.h.s. of eq. ()28dp we observe that the temperature dependent part of the energy flux 
U is the excitation energy flux, Ub = ^^^^^^ Li2 (^^frr) i which is independent of a. For a = 0, 

U = Ub- The microscopic interpretation of this statistics independence is similar to the one 
given in section [2j 

As in Ref. |24j . one can realize a one-to-one correspondence between micro-conflgurations 
of particles that transport the same excitation energies, although they have different exclusion 
statistics. An example of two micro-configurations like this, one of bosons and one of fermions, 
is given in figure [5j Since two such micro-configurations transport the same excitation energies 
and the same particle fluxes, then they transport the same heat fluxes. Moreover, since the 
entropy flux is the flux of the number of micro-conflgurations, the entropy flux should also be 
independent of the statistics |32| . 



4. Conclusions 

I reviewed the thermodynamic equivalence of systems of the same, constant density of 
states (DOS) and the independence of statistics of the heat conductivity of one-dimensional 
(ID) transport channels. I showed that these two classes of apparently different physical 
phenomena have many similarities and practically the formulae that describe the equilibrium 
thermodynamics of constant DOS systems can be directly transcribed for the transport of ID 
channels. 

I also presented the microscopic explanation of these universal features of FES gases. I showed 
that the thermodynamic equivalence is caused by the independence of statistics of the particle 



excitation spectra, whereas the universahty of the heat conductivity of ID channels is due to 
the independence of statistics of particle excitation fluxes. 

In conclusion, the similarity between the mathematical formulation of the equilibrium 
thermodynamics of constant DOS systems and the mathematical formulation of the transport 
properties of ID channels is a reflection of the microscopic similarities between these two types 
of physical phenomena. 
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